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Abstract 

The branching ratio for the rare two-body B s — > fry decay is calculated using the factorization 
assumption. This transition is dominated by the annihilation diagrams and, in principle, prone 
to receiving substantial contributions from new physics. We estimate Br(B s — > pj) 1.6 x 10~ 9 



> 



within the Standard Model and investigate the sensitivity of this decay mode to the effects of two 
new physics scenarios: vector quark model and supersymmetry. Our results indicate that the shift 
j-C in branching ratio is at most around 10% with the addition of vector quarks and is negligibly small 

£N| , in the constrained minimal supersymmetric extension of the Standard Model. 



O 

o 



PACS numbers: 11.30.Pb, 12.15.Mm, 12.60.Jv, 13.20.He 



1 Introduction 



^ ■ With Large Hadron Collider (LHC) becoming operational, search for new physics beyond the Standard 
Model (SM) is entering an exciting new phase. For example, among other particles, a large number of 
B s mesons are expected to be produced at this facility, and as a result, rare processes involving this 
interesting hadronic system will be accessible experimentally. The transition B s —* p-y is of particular 
interest as, unlike B — > pj, it mainly proceeds via the so-called annihilation process. As such, it is 
interesting to study the stability of the SM prediction for the branching ratio of this decay mode to 
the effects of, as yet unknown, exotic particles. 

In this paper, we first calculate the branching ratio for B s — ► decay within the SM using the 
light cone amplitudes and factorization assumption, after identifying the effective operators which 
have the dominant role in this process. To the best of our knowledge, even though our prediction has 
a large uncertainty, this is the first calculation of the branching ratio for this transition. Next, we 
check the sensitivity of this decay mode to possible tree-level flavor changing neutral currents (FCNC) 
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generated by new physics. One particular scenario is the addition of a single down- type vector quark 
to the three quark generations of the SM. Since both the left-handed and the right-handed compo- 
nents of a vector quark are iso-singlets, the weak interactions of this quark can only happen via mixing 
with the ordinary quarks. Consequently, the extended 3x4 quark mixing matrix is not unitary and 
tree-level FCNC is present in this model. In fact, one can show that the non-zero b — ► sZ° in Vector 
Quark Model (VQM) is proportional to U sb = V* s V ub + V* s V cb + V t *V tb which basically measures the 
non-closure of the unitarity triangle within the SM. Our results indicate that, within the acceptable 
range for the above model parameter, the shift from the SM prediction is about 10%. 

We also consider the minimal supersymmetric extension of the Standard Model where every fun- 
damental particle receives a supersymmetric partner with spin differing by 1/2 unit. The particle 
spectrum of supersymmetry therefore consists of three neutral Higgs bosons, two charged Higgs bosons 
and supersymmetric particles such as gluino, squarks, sleptons, neutralinos and charginos. This leads 
to additional contributions to the B s — ► fry branching ratio due to the new particles. We study in 
this paper a scenario with minimal flavor violation and evaluate the numerical implications of such 
contributions. We show that the overall supersymmetric contributions in this decay mode are very 
small. 

We would like to emphasize again that unlike some well known decay channels like the inclusive 
b — > sj, exclusive B — ► K*~f and B s — > <3?7 which have been under thorough investigation in the 
literature, B s — > fry proceeds mainly via annihilation process and in that sense is quite different and 
unique. 

2 Low energy effective Hamiltonian 

The relevant effective Hamiltonian for B s — > fry can be written as 



n -ff = ~m £ v psV Pb [Ci(») o?( M ) + c 2 (aO of M + £ Ci(Ai) Oi(/i)] , 



(i) 




where C's are the Wilson coefficients and the operators are defined as follows: 
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Figure 1: Annihilation diagrams contributing to the B s — > decay. The first row corresponds to A\ 
(left) and Ai (right) and the second row corresponds to ^.3 (left) and Aa (right). 

07 = -^m b s a aL lu PRb a F^ u , 

8 = ^m b s a a^P R T^b G- u . 

The contributions of the current-current operators Of and 0% an( i QCD penguin operators Oi, 
z = 3 ... 6 are shown in Figure [TJ Operators 01, 0\, O7 and Os, on the other hand, have vanishing 
amplitudes to the leading order within the factorization assumption. As a result, the amplitude for 
the B s — ► p7 decay is composed of four generic matrix elements of the following forms: 

A\ = (P( k pi e p)l( k -y, 6y)l*y(l -75)^97m(! ±75)-^ 7i/<w?I-B«(Pb,)) , (3) 

A 2 = (p(fc pi e p )7(fe r ,e 7 )|l7' , (l-76)6?7 I ,^- r - 7m (1 ± 75 )g|S s (p B J) , (4) 

^2 ~~ m g 

•^3 = (p(kp,e p )7(A; 7 ,e 7 )|s7^-— 7 M (1 - 75)^57/^(1 ± 75)<7|#s(pbJ> , (5) 

2> 3 — m s 

A 4 = {p{k p ,e p )-f{k 1 ,e 1 )\s^{l - 75)— — 7^6^(1 ± 7 5 )g|S s (p B J) , (6) 

P4 ~ fnb 

where p\ = k^ + p q , p 2 = k^ + p^, ^3 = p g — kj, and P\ = Pb~ kj, with q = u and <i. 

The contributions from the amplitudes .4i,2,4 are power suppressed as compared to A3 [H[2]. However, 
-4i,2 gi ye a numerically important power correction because of a relative enhancement factor from the 
Wilson coefficients. 

The leading twist light-cone amplitude for a vector meson (V) with flavor content (qq r ) and polarization 
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vector ty can be expressed as p2 H] : 
(V(k v ,ev)\m a q'{z)0) 



z 2 =0 



fvmv 



o 



(7) 



+ 



Jo 



dve mkv - z gy {v \v) 



where the quark and antiquark contents of the vector meson have momenta p q * = vky and pq = vky, 
respectively, with v = 1 — v. The functions f(v) = $y(v), $ v (v), gy (v) and gy (v) are all nor- 



malized, such as / dvf(v) = 1. 
Jo 



The relevant distribution amplitudes for the p meson in B s — * pj are given by [3]: 



(j>y(v) = 6vv 



l + 3aix + 4^X 2 -1) 



(8) 



_L (v) i \ « / 

9v ( v ) = t( 



l + X 2 ) + 4^ X 3 + 04 + 5Cs) (3X 2 -1) 



+ 
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112 °2 + 64 ^i 3 ^ 



(3 - 30 X 2 + 35 X 4 ) 



(9) 



3 ~ 3 ~ 

+- 8 + (2 + In v + In v) + - 5_ (2% + In u - In «) , 



with x = ^ — t> = 2t> — 1. Other coefficients in (|7|), (jHJ) and (fTU|) can be found in Tables [T] and [2j Note 
that in the case of a p-meson, <5 + and 5- vanish. 

The light-cone distribution amplitude for B s meson to leading power can be written as [TJ EJ [6] : 



(0\s(z) a b(0)p\B s ( PBs )) 



z+=zj_=0 



ifB. f 1 



Jo 



where n_ = (1, 0, 0, —1), with the light cone components for any vector v defined as: 

v° ± v 3 



v± 



V2 



(10) 



(11) 



The light spectator carries longitudinal momentum fraction £ = p s+ /pB s +- The normalization condi- 
tions are 



/ d£*B.i(0 = l, /"«*fl.2(0=0 

JO JO 



(12) 



fp (MeV) 


# (MeV) 


of, 


ojj 


a-2 


198 ± 7 


152 ±9 





0.16 ±0.1 


0.17 ±0.1 



Table 1: Decay constants and couplings of distribution amplitudes for p-meson, including SU(3)- 
breaking. The renormalization scale is \j? = 4.8 GeV 2 . 







^3 


T 


u 




Pt 

S4 


0.023 


-1.8 


3.7 


7.5 


0.13 


0.07 


-0.07 



Table 2: Couplings for distribution amplitudes without including SU(3)-breaking. The renormalization 
scale is /j 2 = 4.8 GeV 2 . 



m q 


m p 


m Bs 


3 MeV 


775 MeV 


5367 MeV 



Table 3: Quark and meson masses used in our calculations. 



The contribution of &b s 2 is numerically small [7] and is neglected in the following. We parameterize 
the first negative moment of &b s i by the quantity A# s , such as 

f^ 1 ^, (13) 

Jo £ 

Using the model function from [6], we evaluated the value oi Xb 3 = 580 ± 190 GeV. The numerical 
values for the mass parameters which are used in our calculations are given in Table [3l 



Consequently, from the formulas ([3]) to ([6]), the decay rate can be obtained as: 

2 2 
rvTTlr, — 771., n o 

T(B S - pj) = g B * 3 p G 2 F \V*V tb \ 2 \A\ 2 , (14) 
8 m B 

where \A\ 2 is given as: 

\A\ 2 = |Ci 2 | 2 |^ 1 , 2 | 2 ±|C34|V3,4| 2 + 2Re(C 12 C 3 * 4 ^ li2 ^ i4 ) , (15) 
with Ai t 2 and ^.3,4 defined as the following: 

A lt2 =Ai-A 2 = ± 2f B J p m p {h - h)e KX ^ e K p p Bs k£ ^ , (16) 



A 



3.1 



A? - A 4 



ifB s f p m p 



x 1 



{ Wp ■ fey) (PB. ' £ t) - ( e P • e 7) (PB S • fey)] 



(17) 
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In Eq. (fTUj) . plus and minus signs correspond to V + A and V — A structure, respectively, and Ey in 
Eq. (|17p is the energy of the photon. Eqs. ([7]) and (|10p are used to parameterize the hadronic matrix 
elements. Accordingly, I\ and 1% in Eq. (|16f) are defined as: 



Ii 



o vm 2 B + wrrip — 



Neglecting the small light quark mass, one can show that: 



1 



m 



-5.22 ± 0.33 + (2.47 ± 0.10)X) , 



B a 



0.99 ±0.01 



m 



(18) 
(19) 

(20) 
(21) 



where the logarithmically divergent integral J^dv/v in Eq. (|18|) has been parameterized by X = 
(1 + pe lLp ) ln(mft/A/j) with ~ 0.5 GeV a typical hadronic scale, </> an arbitrary strong-interaction 
phase, and p < 1 [5j. 

C\2 and C34 in Eq. (fl~5|) are combinations of Wilson coefficients appearing in (P): 



C 



12 



3-1 



1 

72 
1 
72 



Ci + 



3 

c 2 



Vus V ub 

v t * s v tb 



+ c 3 + 



c 5 



3 



(22) 
(23) 



The Wilson coefficients are real, however the CKM matrix elements can have imaginary parts which 
are taken into account. The numerical values for C\2 and C34 within the SM are calculated to be 
Cf 2 M = ((0.7 ± 0.1) + i(8.7 ± 0.9)) x 10" 3 and Cf 4 M = ((-1.7 ± 0.1) + i(4.4 ± 0.2)) x 10" 3 using the 
inputs and conventions of [8]. The errors in these coefficients are mainly due to scale dependence and 
determination of the CKM matrix elements. 



Our result for the branching ratio can then be expressed in the following form: 

Br (B s -> fry) = (1.6+™) x lO" 9 ( fB ° ) ' ( Y^t] * . ( 24 ) 

The large uncertainties in the above result are mainly due to the acute sensitivity of our predictions to 
the parametrization of the divergent integral, and therefore this result must be considered as a rough 
estimate for this branching ratio. 

The expected production yield of a decay mode per year at LHCb can be estimated from the fol- 
lowing formula |10j : 

Asignai = 2ff££j£ r x Br(b B s ) x 5r (decay mode) , (25) 
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where a b i, the bb production cross section is expected to be about 500 fib, the integrated luminosity 
for one year of data taking is C l ^ T = 2 fb^ 1 with the assumption of an average luminosity of 2 x 10 32 
cm~ 2 s _1 , and the branching ratio of b quark(antiquark) hadronizing into a B S (B S ) meson is taken to 
be 

Br(b B 8 ) = Br(b -> B s ) = (10.4 ± 1.4) x 10" 2 . (26) 

Based on Eq. (|25l) and our estimated branching ratio in Eq. (|24l) . one would expect about 340 
B S (B S ) —* p7 per year. The signature of this decay mode is two pions with invariant mass around p 
resonance plus a photon. The key in successful measurement of this decay channel could be the high 
energy resolution of the calorimeter so to make it possible to make the required background cut which 
is expected to be mainly due to B —* p7 process. The smallness of the likelihood of this process as 
expected by the SM could make it quite sensitive to the new physics, especially those which allow 
FCNC at the tree level. We examine one such beyond the SM scenarios as well as supersymmetric 
contributions in the next section. 



3 Contributions from Vector Quark Model and Supersymmetry 

In the vector-quark model (VQM) |llj . in which SM is augmented by an additional isosinglet down 
type quark the 3x4 quark mixing matrix is not unitary, which leads to non-vanishing FCNC at the 
tree level. For example, the elementary b — > sZ° vertex is proportional to the non-unitary parameter 
U sb = (y^V) sb at the leading order which directly measures the non-closure of a unitary triangle and 
as such, it is quite desirable to find observables which are most sensitive to this parameter. In fact, 
one can show that the factor C±2, as defined in (I22p . receives additional contributions 1121: 



SO 



12 



1 

V2 



VQM up 



where 



a 



VQM 



a 



VQM 



u 



sb 



VtbV t 



n VQMup\ . 1 / n VQM down 
U 5 J+ol°3 



-{I q w - Q q sm 2 9) 



a 



VQM down\ 



(27) 



(28) 



ts 



U 



sb 



v tb v t * s 
u sb 



n 2 2 

sin 6 

2 3 

1 1 , 
1 — sin' 

2 3 



= 0.35 . . . q = up 

9 = -0.42 . . . q = down 



VtbVts 



-Qqsm 9 



(29) 



U 



ah 



v tb v t * s 
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-— sin 
3 

■ sin 2 9 



= —0.15 . . . q = up 
0.08 . . . q = down 



One can use the experimental measurement of the dileptonic B — > X s £ + £~ decay by BELLE [13] 
and BABAR p3] to constrain the acceptable domain for the magnitude and phase of U sb . It turns 
out that this parameter can attain a maximum value of 1.1 x 10~ 3 if its phase is close to zero [15] . 



The resulting shift in the C\2 coefficient in this model is at most 8C 



VQM 
12 



3.1 x 10 , leading to 



a branching fraction of about 1.8 x 10 . Therefore the SM branching ratio in (124H is increased by 
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around 10% with the addition of vector quarks. 



We also consider super symmetric contributions to the B s — * p-y branching ratio. In particular di- 
agrams involving charged Higgs boson, charginos and squarks can modify the Wilson coefficients 
and consequently the branching ratio. To study the numerical implications of such contributions we 
restricted ourselves to the constrained minimal supersymmetric Standard Model (CMSSM) - often 
called mSUGRA - with the number of free parameters reduced to five assuming the unification at 
GUT scale. In this scenario, the chargino-squark contributions appearing in the loops are completely 
negligible. The charged Higgs boson on the other hand can replace the W boson but its contribution 
is still smaller than the SM contributions. In fact, within the acceptable range for the 5 parameters 
of mSUGRA, the upper bound for ^(7 mSU GRA j g sma n er than the SM contribution by 3 orders of 
magnitude and therefore, one can safely assume that the prediction in (I24|) is not affected by the 
constrained SUSY particles. 

4 Conclusion 

The rare decay B s — * frj could be an interesting process to look for at the LHC. It may be an 
experimental challenge to measure this decay mode due to the fact that the final state particles are both 
neutral. However, knowing the branching ratio helps to devise a strategy for eventual experimental 
measurement. 
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